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Abstract 



In a time dependent background like de Sitter space, Feynman- Dyson perturbation the- 
ory breaks down due to infra-red divergences. We investigate an interacting scalar field theory 
in Schwinger-Keldysh formalism. We derive a Boltzmann equation from a Schwinger- Dyson 
equation inside the cosmological horizon. Our solution shows that the particle production 
is compensated by the reduction of the on-shell states due to unitarity. Although the de- 
grees of freedom inside the horizon leads to a small and diminishing screening effect of the 
cosmological constant, there is a growing screening effect from those outside the horizon. 
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1 Introduction 



Investigating field theory in de Sitter space may illuminates deep mysteries surrounding 
inflation in the early universe and dark energy of the present universe. The past and current 
exponential expansions of the universe are likely to be driven by the effective cosmological 
constants of the order of GUT and neutrino mass scales respectively. We are perplexed by 
the huge disparity of the relevant energy scales. Of course we do not understand why they 
are small in comparison to the Planck scale in the first place. Phenomenologically it appears 
that the cosmological constant has evolved with time. Although we may parametrize it by 
a scalar field with a suitable potential, its microscopic understanding is totally lacking. 

Slow roll inflation models possess approximate conformal invariance and the conformal sym- 
metry plays an important role to understand the magnitude of the correlators [1] [2] [3] and 
possible dS/CFT correspondence [1][5][6] . In string theory, there seems to be no stable de 
Sitter vacuum as we need to consider brane-antibrane systems to realize it. 

Since our understanding is so sparse, we wonder if we are entering a completely new territory. 
It might very well be the case since the standard Feynman-Dyson perturbation theory breaks 
down in a time dependent background like de Sitter space. Feynman-Dyson formalism is 
the backbone not only in relativistic field theory but also in classical statistical mechanics 
and critical phenomena. In this sense our expertise might be confined in equilibrium physics 
while our problem belongs to non-equilibrium physics. 

In fact we need to use Schwinger-Keldysh formalism to investigate field theory in a time 
dependent background like de Sitter space. We can derive a Boltzmann equation in this 
formalism which is a standard tool to investigate non-equilibrium physics [7]. In such a 
setting, it is in principle possible that the effective cosmological constant changes with time. 
In other words the dynamics may explain deep mysteries of this century if we can demonstrate 
that the cosmological constant decreases with time in an interacting field theory. Although 
there exist several proposals along this line of thoughts in the literature, our understanding 
is still in a preliminary stage [7] |8] [9] [lO] [11] . 

There is a long history of studying Boltzmann equations in Schwinger-Keldysh formalism 
[12] starting from Kadanoff-Baym [13] [S] [IS] . In this paper we derive a Boltzmann equation 
in de Sitter space from a Schwinger- Dyson equation. This problem has been studied to the 
leading order of the derivative expansion of the Moyal product in the Wigner representation 



1 



[20] . However only the energy conserving process has been identified in such a hmit. We go 
beyond the leading order of the expansion to investigate the particle production effects due 
to energy non- conservation in de Sitter space. We also investigate the energy-momentum 
tensor of an interacting scalar field theory to estimate the effective cosmological constant. 

The organization of this paper is as follows. In section 2, we introduce a scalar field theory 
in de Sitter space. In section 3, we recall Schwinger-Keldysh formalism. In section 4, 
we determine the full propagator inside the cosmological horizon to the leading order in 
perturbation theory. In section 5, we estimate the effective cosmological constant from the 
energy-momentum tensor of a scalar field. We conclude in section 6 with discussions. 



2 Scalar field theory in de Sitter space 

In the Poincare coordinate, the metric in de Sitter(dS) space is 

ds^ = -de + e2^*dx2 

-rfr2 + dx2 (2.1) 



where the dimension of dS space is taken as D = 4. if is the Hubble constant and the 
conformal time r is related to the cosmic time t as r = —jje~^^. It assumes the value in the 
range — oo < r < and it increases with cosmic evolution. 

In this paper, we consider a massless scalar field (p which is minimally coupled to the dS 
background. The quadratic action for the matter field is 

Smatter =\ J V^d^X [-g^'' d^ifd^ifi]. (2.2) 

We find it convenient to redefine the scalar field as follows (f — )■ Hnp. We can simply scale 
it back to find the original scalar field. In terms of the rescaled field, the quadratic action 
becomes 

Smatter = ^ j d^X ^ (^-9^ + + (2.3) 

The positive frequency solution of the equation of motion with respect to the above action 
is 

0p(a;) = 4=(l-^-)e-*'^^*P^ (2.4) 
where p = | p | . We expand the scalar field as 

/d^p 
(ap0p(a;) + aj,0;(x)) . (2.5) 



We consider the Bunch-Davies vacuum |0) which is annihilated by all the destruction oper- 
ators Vap|0) = 0. The propagator in such a vacuum is 

(2.6) 



(1-i — )(l + z — ) e- 



(27r)3 2p pn pT2 

If we blindly apply Feynman rules to investigate the effects of the interaction, the integrations 
over time give rise to infra-red(IR) divergences at the infinite future [TB]. For example, with 
X{p^ interaction we find 

r dr, -1^ X {iHT)iHn){T^ixMx^))Y ~ [' dn (2.7) 

In the next section, we recall the method to investigate the effects of the interaction which 
does not suffer from these IR divergences. 



3 The Schwinger-Keldysh formalism 

Let us represent the vacuum at t — ?■ — oo as |m), and t — +oo as \out). In the Feynman- 
Dyson formalism, the vacuum expectation value(vev) is essentially given by the transition 
amplitude between \in) and \out) 

(Onit)) = {out\U{+oo,t)Oi{t)U{t,-oo)\in) 

= {out\T[U{+oo,-oo)Oj{t)]\tn), ^^'^^ 

where Oh and Oj denote the operators in the Heisenberg and the interaction pictures re- 
spectively. U(ti,t2) is the time translation operator in the interaction picture 



U{h,t2) = T 



exp|-i^ dt Hiit)"^ 



(3.2) 



It is because \in) is equal to \out) up to a phase due to the time translation invariance. On 
the other hand, there is no time translation symmetry in de Sitter space, and so we can't 
prefix \out). In this case, we can evaluate the vev only with respect to |m) 

{Onit)) = {in\U{-oo,t)Oi{t)U{t, ~oo)\in) 

= {tn\Tc[UOi{t)]\m). ^^'^^ 

Here Tq denotes the operator ordering specified by the following path 
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+00 



dt= dt+- I dt-. (3.4) 

' J — oc J —oo 

Because there are two time indices (+, — ), the propagator has 4 components 

^ " - [G-^{Xi,X2) G--{x^, X2) J 



^ f {T^{xiMx2)) {(p{x2)(f{xi)) 



(3.5) 



Here T denotes the antitime-ordering. 

In order to investigate the effects of the interaction, we consider the foUowing Schwinger- 
Dyson equation for the two point function 

G{Xi,X2) = Go{Xi,X2) 

+ J ^/^d'^X3^/^d'^X4 Go{Xi,X3) ^gg-j 

X t{x3,X^) G{X4,X2), 

where Go is the free propagator, G is the full propagator, and E is the particle's self energy. 
Especially, we focus on the ( — h) component of the propagator 

G''"^{xi,X2) = G^'^{xi,X2) (3.7) 
+ J \/-gz'i''xz\/-gAd''x^ G'o+(xi,a;3)E++(a;3,X4)G'++(x4,X2) 

- ^ V-93d^X3y/^gld'^X4 G^'^{xi,X3)E^~{x3,Xi)G~^{x4,X2) 

-J V-g^d'^Xs^Z-Qid^Xi Go~(xi,X3)E"+(x3,X4)G'++(x4,X2) 

+ J ^/^d'^X3^/^d'^Xi Go~{Xi,X3)T,~~{x3,X4)G~^{Xi,X2) 
= Gq'^{xi,X2) 

+ J \f^zd^X3sJ-g4d^X4 Gq{Xi,X3)T,^{x3,X4^)G'^{x4,X2) 

+ J -g-id^X3sJ -g^d^x^ G^(xi,X3)E~+(x3,X4)G^(x4,X2) 

+ y -gzd''x3sj -gd^'x^ Gq'^{xx,X3)YA[x3,X4)G^[x4,X2). 
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Here we have introduced the retarded and the advanced propagators as follows 

G^txi, X2) = e{t, - t^) [G'-+(xi, 0:2) - G'+-(a;i, X2)], 
G^(xi,X2) = -e{t2 - ti)[G'-+(xi,X2) - G'+-(xi,X2)]. 

We observe that the integrations over time are bounded by ti or t2 because of the causality. 
In the same way, the following identity also holds 

+ j V-gld'^X'iV-Fid^Xi G'^(Xi,X3)S^(x3,X4)G'o+(x4,X2) 
+ j V^d'^Xsy/^gld'^Xi G'^(xi,X3)S" + (x3,X4)Go (X4,X2) 

+ j V-g3d^X3y/-g4d^Xi G"+(xi,X3)S^(x3,X4)Go (3^4,3:2). 

In this formalism, the integrations over time are manifestly finite due to the causality. This 
formalism is called Schwinger-Keldysh formalism. In order to understand the effects of the 
interaction, we derive a Boltzmann equation on the dS background from a Schwinger-Dyson 
equation in the next section. 



(3.9) 



4 Boltzmann equations from Schwinger-Dyson equa- 
tions 

In a time dependent background, we need to consider excited states in general. For such 
a state, the expectation value of the number operator (a'^a) is non-vanishing. We introduce 
a distribution function / for scalar particles as follows 

(a;aq) ^/(p)x(27r)35(3)(p-q). (4.1) 

One of our main objectives in this paper is to understand the time dependence of the dis- 
tribution function / due to the interaction. We utilize a Boltzmann equation for this pur- 
pose. Boltzmann equations govern the time evolution of the distribution functions. They 
are widely used to study non-equilibrium physics. In fact there is a long history of the 
microscopic derivation of Boltzmann equations in non-equilibrium physics using Schwinger- 
Keldysh formalism |13] [13] [IS] . In this paper, we systematically investigate the propagator 
in dS space from a Schwinger-Dyson equation. 
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(4.2) 



We assume that the full propagator in de Sitter space has the following form 

The propagator depends on the average and the relative time: 

Tc = — ^ — , r = Ti - T2. (4.3) 

It consists of the on-shell part and the off-shell part. In the on-shell part, we have introduced 
the wave function renormalization factor Z{p, t^. The off-shell part depends on the spectral 
function F±{e, p, r^). We assume that /, Z, F± evolve with the average time Tc- We investigate 
the propagator in the region: 

kcl > |re| > 1/p. (4.4) 

The second assumption implies that we investigate the propagator well inside the cosmolog- 
ical horizon. 

From (13. 7p and (13. 9p . we can derive the following identity 

Gq^\iG ^{Xi,X2) — Gq^\2G '^{Xi,X2) 
= + V-gi j V-g3d^X3 E^(X1,X3)G'"+(X3,X2) 

+ J V^d'^X3 J:'^{xi,X3)G^{x3,X2) 

- yf^2 \ V^d'^X3 G^(X1,X3)S"+(X3,X2) 



- \/^92 J V^gsd X3 G +(xi,X3)S (X3,X2). 

By putting the expression for the full propagator (14.21) into the left-hand side of the Schwinger- 
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Dyson equation (14. 5p . we obtain 

Gq^\iG ^ {Xi, X2) — Gq^\2G ^(Xi,X2) 



+ 



(27r)3 l^dTc pdrdr, 
decPp 



e>0 



(2vr) 



-ier+ip-x 



Here we recall the following definitions 

G-,'^t{dl-dl-^\ 

G-^\G\xr.X2) = b'^''\x^-X2), 
G-^%G^{x^,X2) = ^^""Kxr - X2). 
In (14.61) we have shown the leading terms in the power series expansion of l/pTc 



(4.6) 



(4.7) 



The right-hand side of Eg. (14.51) corresponds to the collision term G\f\. In this paper, we 
investigate the effects of the interaction in Ay?'^ theory at the one loop level. We subsequently 
find that this theory captures the essential features of more generic field theories such as g'^'^ 
theory. The self-energy is 




S*-'(X3,X4) 



-zA)^ 



G*^(X3,X4)G"^(X3,X4), i,j = +, 



(4.8) 



To the leading order in perturbation theory, we can approximate that /(p, Tc 
Z{Pj'^c) = 1, F±{£,P,^c) = in the collision term. We also expand the collision term 
by the power series in l/|pTc| type factors which can be justified well inside the cosmological 
horizon. It is a kind of the derivative expansion of the Moyal product in the Wigner repre- 
sentation. We indeed find the particle production effects due to the non-conservation of the 
energy in this expansion. 

In this investigation, we need to perform the following integrations at the interaction vertices. 

1 



3 

-00 '3 



i{e±p)T3 



neN, i = l,2, 



(4.9) 
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where e — ±pi ±p2- We evaluate these integrations in the assumption \{e±p)Ti\ ^ 1 . For 
our purpose, it suffices to evaluate them to the next leading order 



n 1 



-n 



(4.10) 



By using these approximations, we derive a Boltzmann equation in de Sitter space. In what 
follows, we investigate the collision terms and their properties in detail. 



We henceforth suppress the following integration factor in the propagator 



/ 



(4.ii; 



In other words we work in the momentum space by performing the Fourier transformation 
with respect to the spacial coordinate x. 
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4.1 The structure of the colhsion term 



From the Schwinger-Dyson equation (14.51) . we observe that the colhsion term has the 
on- shell part and the off- sell part. Firstly, the on-shell part comes from the following contri- 
butions 



Con[f] = + y/^9i j V-gsd'^xs J:^{xi,X3)G +(0:3, 2:2) 



oc e 



We evaluate the on-shell part to the leading non-trivial order (9(l/r^) as 



Conlf] 

= -(l + /(p))e 



-ipT 



^ de 1 1 Jf 



X / ' dp,{il + fip,m + fi6-p,))-fip,)fi6-p,)} 

r de . 1 1 if 1 1 -K 

X 2^^^7^ ^ 7+^'^ ^\6-py ~ {e + py'^^ 

POO 



2 



+ f{p) e+^^" 



■X 



1 ^ 1 

271^ e — p e + p' t'^ ' "^{e — py i^ + P^^'^c 



-IT , 1 

+ 



} 



2 



+ 2 



"Pder. I 1 .-if , 1 1 .-1 



:{( 



27r ^ £ — p e + p 



} 



X + /(pi))/(pi - £) - /(pi)(i + /(pi - £))} 

See Appendix A for the details of the calculation. 



(4.12) 



(4.13) 
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Secondly, the off-shell part originates from the following contribution 



(4.14) 



rpi+p 

dpi / dp2 e-*(±fi±P2)^ 
"'Ipi-pI 



"OO 

OC / UPl 



The off-shell part is also calculated to 0(l/r^) as 

Cofflf] 



_l X 

IGvrp^iJ^ 
de 



(4.15) 



27T 



^ ^Zl f ' dpi (l + /(pO)(l + /(e-pi)) 



+ 2 



2^ 



c^Pi {l + fiPi))fiPi-e) 



IGnp'^H'^ 
de 



-1 f 2 



27r 



+ 2 



2^ 



dpi f{pi)f{e-Pi] 



dpi f{pi){l + f{pi-e)) 



We note that the both on-shell fl4.13p and off-shell fl4.15p collision terms have infra-red 
divergences aX e = p. There is a standard procedure to deal with this problem in massless 
field theory and we find that it also works here. First of all, we need to recall that any 
experiment has a finite energy resolution Ae:. So we need to add the on-shell and off-shell 
collision terms within the energy resolution Ae:. We first divide the integration range of 
Coff[f] as follows 



oo /"oo />p+Ae 

+ / 



Jo Jp~Ae 



(4.16) 



We then redefine the on-shell term C'^^[f] and the off-shell term C^jj[/] by transferring the 
contribution of Co//[/] within the energy resolution p — Ae < e < p + Ae to Con[f]- The 
explicit expressions are shown in Appendix A. 

When f{p) = 0, we find that infra-red divergences cancel out in this procedure. In the 
next subsection, we investigate the case when / is a thermal distribution. For a generic 
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distribution, the cancellation does not take place and we seem to face linear IR divergences. 
However there is no real infra-red divergence in our problem since the time integration range 
in f l4.10p is bounded by r^. We thus argue that the linear divergence should be cut-off at 
~ 1/\t^\. 

Before investigating the thermal distribution case, we point out the difference between 
Minkowski space and dS space with respect to the collision term. In Minkowski space, 
the collision term does not have the off-shell term due to the time translation symmetry 

CofAf] cx 1^ 2ir5{e - p)e^-' = e^^^' =^ C'^^fif] = 0. (4.17) 

On the other hand, as we observe in fl4.15p . the collision term in dS space has the off-shell 
term due to the absence of the time translation symmetry. This is why we have introduced 
the spectral function F±{e,p,Tc) in the full propagator fl4.2l) . 



4.2 Thermal distribution case 



We focus on the case that the initial distribution function is thermal in this subsection 

1 



fip) 



(4.18) 



where /3 is an inverse temperature. In Minkowski space the thermal distribution is obtained 
as the solution of the Boltzmann equation. On the other hand, we find that the collision 
term in dS space is non-vanishing even for the thermal distribution. 



The off-shell collision term can be evaluated as follows 

C'offif] 



de 
2^ 



(l + /(^))e- 



;i + G(e,p,/3)) 



p+Ae 

Jo 27r {e-pf {e + pf rl 

X 



Vo'KpH'^ 

°° de 1 



p+Ae 



2tt 



(e-py {e + p^ rl 



1 ^^(l + G'(£,p,/3)) 



^P-'^^ de , , 







27T 



^^G{e,p,f3) 



(4.19) 
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where 

G(.,p,«.Alog(i^^|. (4.20) 



We note that the above expression is of the following form 

P dF 

C'off[f] = / -((l + /(e))A(5,P,r,)e-^-/(eK(5,p,r,)e-^). (4.21) 
It is consistent with our ansatz for the full propagator f l4.2p . 
Finally the on-shell collision term is evaluated as follows 

C'on[f] 

(l + /(p))e-*^"x (4.22) 



d£ ^ ^ 1 ^ 1^ _ 1^ -1 

r-^'de.. 1 1 .if . 1 1-1.. 

+ 1 ^^^7^p + 7T-p^-! + ^(F^ - (FT]^)^^^^^'^' . 
^{(^ + ^)^ + (r^2-r^2)^}(^ + G{e,p,p)) 

r-^'der, 1 1 .-if . 1 1 .-1.^. 

The details of its derivation can be found in Appendix A. 

Here we have cut-off the IR log divergences when \e—p\ ~ 1/|tc| because our time integration 
f l4.10p does not diverge even when e = p. From the on-shell collision term f l4.22p . we observe 
that it is necessary to introduce the wave function renormalization factor Z{p,Tc). In the 
last line, we find that the remaining logarithmic IR contribution leads to the modification 
of the thermal distribution function 6f{p,Tc). 

So far, we have focused on the IR singularities due to the interaction. Of course, there are 
also the ultra-violet (UV) divergences in the collision term. The off-shell part f l4.15p does 
not have the UV divergences because of the exponentially oscillating factor. We also assume 
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that a generic distribution function vanishes exponentially at the UV region like the Bose 
distribution 

f(P^) - ^ 0- (4-23) 

From these facts, the UV divergences in the collision term is estimated as follows 

C[f]uv = C'^if] 



+ 16^^^^^' 2^^7^^7T^^^ (4-24) 

lOTT^ -n Tc <? 2p 

Since the integral is logarithmically divergent, we need to introduce a UV cut-off. We argue 
that we need to cut-off the integral at a fixed physical energy scale Auv- As the physical 
energy is eH\T\, this prescription leads to a time dependent UV cut-off Ai/y/Hlrl = Kjjye^^ 
in the above expression. We believe that this is a physically sensible prescription which is 
consistent with general covariance. In this prescription, the degrees of freedom inside the 
cosmological horizon remain the same with respect to time. The IR cut-off is provided by our 
energy resolution Ae in fl4.24p as the IR singularity is canceled by the off-shell contribution. 
The final expression logarithmically depends on the virtuality = {p + Ae)^ — . 

This UV divergence is renormalized by introducing a mass counter term in the action which 
leads to the following collision term 

C[f]s^^ = + ^^,Sm' X (1 + fip))^e-''" 

-^J^^^'x/(^^)^e^^''' (4.25) 

.2 

om =— -log , 

lovr^ /i 

where /i is the renormalization scale. After the renormalization, we obtain the following 
effective mass 



m 



eff IQ^2 



/ a 1 r°° 1 1 /I _ p-/3(£+P)/2\ \ 



including the finite temperature correction. In the zero temperature limit, it agrees with the 
renormalized mass in the flat space. 
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The IR logarithm in the colhsion term fl4.22p leads to the change of the distribution function 
as we solve the Boltzmann equation 

A2 1 



1 



^\og\AeT,\np) 



log lAerr, 



/3 



e/3p _ 1 e/3p - 1 ' 



The wave function renormalization factor is determined as 

X 



327rpiJ2 

°° de 



+ ■ 



,0 27r (e-p)^ (e + p)^ 
The off-shell part of the propagator is determined in terms as 



= + 



327rpif2 



(l + /(^))x 
1 



[e-pf [e^pf tI 



-rd{p-e){- 



1 



F_(e,p,rc) 
= + 



?,2TipH'^ 
' 0{e-p){ 



1 



1 



1 



G{e,p,(3) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



(e-py {e + pY rl 

We observe that the on-shell weight represented by the wave function renormalization factor 
Z is reduced from the unity in a consistent way with the off-shell spectral weight. In this 
sense unitarity is respected by the interaction. 



We have thus determined the full propagator inside the cosmological horizon to the leading 
order of the perturbation theory. We have found that the full propagator which is character- 
ized by f lOTj) . dSH]), (SSnD, fOO]) depends on r^. At first sight, it appears to change with 
cosmic evolution. More and more off-shell states are created with a lapse of time as on-shell 
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states are correspondingly reduced. However we may represent (I4.27p . fl4.28p . fl4.29p . fl4.30p 
by the physical quantities, 

X = -^, P = H\t\p, AE = H\T\Ae, T = H\t\1, M = H\t\^i. (4.31) 
n\T\ p 

In terms of the physical quantities, the full propagator of the original scalar field at the equal 
time f = is 

.2 



G ^{Xi,X2) 



(l + 2(/ + 5/))(l + (5Z) 1 + (1- 



(27r)32P 



m 



2^2 ^ p2 ; 



X 



m 



eff 



Q 



2T 



log-^ / dE\ 

P I ^ 



+ 



E-P E+P 



}logf — 



1 _ ^-{E+P)/2T\ 
e-\E-P\/2T J 



5f 



X' df{P,T)^ AE 



647r2p dP 



log 



H 



6Z = -^x 
327rP 

1 . 1 



1 



dE 



. 2tc AE 2P Jp+AE 

271 \E-py {E + py P 

^-^^ dE 



c^rp I 1 _E±P 

21 ^ I 1 — e 2T 
log 



1 — e 2T 



+ -TT- 



2TT\E-py {E + py P 



2r, 







E+P \ 






— e 


IT 


] 






P~E 




— e~ 


2T J 





(4.32) 

(4.33) 
(4.34) 



(4.35) 



327rP 



:i + fiE,T))x 



1 1 2T 



E+P 

1 — e 2T 
1 — e 2T 



+ ^(P-^)( 



{E-py {E + pf p 



2r, /l-e-^ 

— log I 3p3^ 

1 — e 2T 



(4.36) 



327rP 



f{E,T)x 



+ ^(p-^)( 



1 



{E-Pf {E + Pf P 



2r, 

-FT log 







E+P \ 






— e 


2T 








P~E 


] 




— e~ 


2T , 





(4.37) 



We find that the explicit Tc dependence disappears in these expressions. If we focus on the 
physics at the fixed physical energy scale E, it remains the same with cosmic evolution. It is 
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a very sensible conclusion as we do not expect physics such as particle mass to change with 
cosmic evolution. For a fixed e, the physical energy E decreases with time evolution. So the 
cosmic evolution is identical to the evolution under the renormalization group. We recall 
here that the radial coordinate in AdS space corresponds to the energy scale in AdS/CFT 
correspondence. Since the radial coordinate in AdS space is related to the time coordinate 
in de Sitter space by analytic continuation, de Sitter space seems to be related to AdS space 
in this respect. The only physical time dependence appears through the temperature T as 
it cools down linearly with Tc for a fixed /3. 

We do find a non-trivial modification of the distribution function from the Bose distribution 
due to a large IR effect. The effect of the interaction on the distribution function f l4.27p is 
such that it reduces the particle density in comparison to the Bose distribution. This effect 
can be understood as follows. A single particle can turn into two particles due to the cubic 
interaction. So such off-shell two particle states are created while the on-shell state weight 
is reduced by the same amount due to unitarity. The off-shell states cost more energy and 
so are less numerous due to the Bose distribution function. The net effect is the further 
reduction of the particle density. 

In this section, we have investigated the effects of the interaction on the propagator well inside 
the cosmological horizon. The spectral weight of the off-shell states increases with time while 
the weight of the on-shell states decreases due to the interaction. The modification of the 
Bose distribution is analogous to QCD where the logarithmic divergence requires the scale 
dependent modification of the parton distribution function. In term of the physical energy 
and momentum variables, explicit time dependence disappears and the time evolution may 
be identified with the renormalization group evolution. So we find that the effects of the 
interaction in de Sitter space parallel to those in flat space. As it is explained in Appendix B, 
these features are also shared by gip^ theory. We thus expect they are the universal features 
of the interacting fleld theories in de Sitter space. 

Nevertheless we should keep in mind that we have investigated the propagator near flat space 
and the expansion in terms of 1/pTc breaks down near the cosmological horizon. To fully 
understand the behavior of the two point function in de Sitter space, we have to extend our 
work to the region \pTc\ ~ 1 and \pTc\ ^ 1. 

Before concluding this section, we briefly investigate the non-thermal distribution case. The 
modification of the distribution function 6f{p, Tc) by the cubic interaction is roughly of the 
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following magnitude 

dfip,Ti,T2) 



1 



^"^C lp+\l/Tc\ ' 

(4.38) 



A2 1 



P H'^Tc 

So it is 0(1/p|tc|) in a generic case instead of 0(l/(pTc)^) for the thermal case. While it is 
much larger than the change of the thermal distribution when p\tc\ ^ 1, it becomes only 
important near the cosmological horizon. Although the thermalization may take place when 
the couphng is strong enough A > if, it could only occur near the cosmological horizon. At 
the higher loop level, the thermahzation could also take place through the effective n point 
couplings. We find this is a very interesting problem which requires further investigations. 



5 Effective cosmological constant 

In this section, we investigate the energy-momentum tensor of a scalar field. Since it 
appears on the right-hand side of the Einstein equation, it provides us an important clue to 
understand the back-reaction to de Sitter space. The cosmological constant is renormalized 
order by order in perturbation theory as well as other coupling constants such as the Newton's 
constant. If these renormalization effects are time independent, we can assume that they 
are canceled by the counter terms. We are thus most interested in whether there is a time 
dependent back- reaction to the cosmological constant. 

The energy-momentum tensor for a minimally coupled scalar field is 

2 SSjjiatter 



(5.1) 



The Einstein equation is 

Rf,u -^g^^uR + ■^g^l,. = i^T^u, k^SttG, (5.2) 

where A is the cosmological constant and G is the Newton's constant. The effective cosmo- 
logical constant Agff may be estimated as 

Ae// ^ A - It/. (5.3) 
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The vev of T^'^ may be evaluated as follows 

4A 

Pi 

(5.4) 



where we have used the equation of motion. 

Let us investigate the expectation value of the energy-momentum tensor for a free massless 
field case first. We may formally evaluate the energy- momentum tensor as 

^ P I 2 tt2 2 , 



Although they are quartically divergent, the UV contribution to the renormalized energy- 
momentum tensor is determined from the conformal anomaly T^^ |T7j [18] 

1 29 

We may assume that its contribution to the effective cosmological constant through fl5.3p is 
canceled by a counter term in the bare cosmological constant. From the Bose distribution 
function /, we find the following contribution 



Too = ii^T'' + -H'T')—— = E 



30 24 {Hrf {Ht 



i2' 



2 1 1 1 

= (— T''-—H^T^)-^ = V-^, i = 1,2,3. (5.7) 

Since the physical temperature T ~ |r|, the energy £ and pressure V decreases with time 
and the cosmological constant becomes more and more dominant once T < H. 



5.1 Contribution from inside the cosmological horizon 

Firstly, we evaluate the contribution to the energy-momentum tensor from the degrees 
of freedom well inside the cosmological horizon |pr| ^ 1. From the two point function, we 
can estimate the trace of the energy- momentum tensor as in (15.41) . We substitute the full 
propagator specified by flCTl) . f lCTD . fICTD . f lOOD into -| (^^) 
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Of course, we can obtain the same result without solving the Boltzmann equation. We have 
checked that the following expression in Schwinger-Keldysh perturbation theory reproduces 
the identical result 

= (ml (^)7_1^ fiw^^ I ^'^^^'^^ mx),^%x,)],^^x,)] |m). 

(5.9) 

We also need to estimate the contribution from the three point function in f l5.4p . It turns out 
that the three point function does not contribute to the effective cosmological constant as it 
is O(A^T^). See Appendix C for the details. So the contribution to the effective cosmological 
constant from the degrees of freedom inside the cosmological horizon is 



— K- 



Although we find that it reduces the original cosmological constant, it vanishes as the universe 
cools down with time. It therefore does not lead to a decreasing cosmological constant. In 
order to find such an effect, we may need to consider the contributions to the cosmological 
constant from the degrees of freedom outside the cosmological horizon as follows. 



5.2 Contribution from outside the cosmological horizon 

In this subsection, we calculate the contribution to the energy-momentum tensor from 
the degrees of freedom outside the cosmological horizon: |pr| ^ 1. Inside the cosmological 
horizon, the degrees of freedom are constant because we adopt a fixed physical UV cut-off 
A and IR cut-off H corresponding to the following energy integration range 

.A/H\r\ 



i/\t\ 27r 



(5.11) 



On the other hand the degrees of freedom outside the cosmological horizon increase as time 
goes on. They correspond to the following energy integration range 

de 

^, (5.12) 

where Eq denotes an IR cut-off. Simply put, the degrees of freedom outside the cosmological 
horizon increase as more and more degrees of freedom go out of the cosmological horizon 
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with cosmic evolution in de Sitter space. Since every degree of freedom can contribute to 
the cosmo logical constant, they may give rise to an accumulating effect. In fact we find 
that (T^^') receives a time dependent contribution from the degrees of freedom outside the 
cosmological horizon. 

To investigate this effect through the trace of the energy- momentum tensor (15. 4p . we need 
to calculate the expectation value of the two point and the three point functions. As it is 
explained in Appendix C, (T^^) is estimated for the / = follows 

(T/) = -\ (V^^) - ^(V^^) = ^ ^(-log(-.or))^ > 0. (5.13) 

We thus find that the contribution from the degrees of freedom outside the cosmological 
horizon gives rise to a very interesting effect. It screens the cosmological constant as the 
universe evolves 

Ae// = A-^^(-log(-eor)f. (5.14) 

In this paper we have investigated quantum effects which diminish the cosmological con- 
stant. Of course a scalar field can classically roll down a potential to lower the effective 
cosmological constant. Therefore it is interesting to ask whether quantum screening effects 
exist in classically stable theory such as g^p^ theory. In this case, (T^^) is given by a two 
point function after using the equation of motion 

{T,n = -\v'{^')- (5.15) 
The expectation value of the two point function is 

<^'^ = 2e.'33''L6 (-M-^or))^ (5.16) 
The effective cosmological constant is estimated as 

Ae// = A - -^^(-log(-eor))^ (5.17) 

It turns out that the cosmological constant is also screened in gip^ theory in an analogous 
way . Our findings in this subsection are in agreement with |19] . 



6 Conclusions 



In this paper, we have investigated an interacting scalar field theory in de Sitter space. 
As Feynman-Dyson perturbation theory breaks down, we need to employ Schwinger-Keldysh 
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formalism. Our problem therefore belong to non-equilibrium physics. It may illuminates 
great mysteries of our time: inflation in the early universe and dark energy today. Phe- 
nomenologically it appears that cosmological constant is not a constant but it changes with 
cosmic evolution. Since the vacuum changes with time in de Sitter space, it is logically 
possible that the effective cosmological constant also changes with time. It is because the 
expectation value of the matter energy-momentum tensor contributes to the effective cos- 
mological constant through the Einstein's equation of motion. We are thus most interested 
in a possible time dependence of the effective cosmological constant. 

We have first investigated the time dependence of the propagator well inside the cosmological 
horizon. We have derived a Boltzmann equation from a Schwinger-Dyson equation. We have 
found that the total integral of the spectral weight remains to be unity as the particle creation 
effects are accompanied by the reduction of the on-shell states. The leading IR contributions 
cancel between the real and virtual processes. This fact suffices for the complete cancellation 
of the infra-red divergences at zero temperature as pointed out in section 4. However it 
is not so at finite temperature. The remaining IR contribution leads to the modification 
of the particle distribution function as f l4.34p . This interesting effect vanishes in the zero 
temperature limit as the Bose distribution function itself vanishes. 

Although the propagator is time dependent, explicit time dependence disappears when it is 
expressed by physical quantities. In another words, time evolution may be identified with 
the renormalization group evolution in an interacting field theory in de Sitter space. We have 
thus found that a field theory in de Sitter space shares many common properties with those 
in Minkowski space. We believe this is due to unitarity of the theory. We also believe it is 
due to the fact that the degrees of freedom inside the cosmological horizon does not change 
under cosmic evolution. In any field theory we need to adopt a fixed physical UV cut-off 
A while IR cut-off is provided by the Hubble constant. From this reason, we find that the 
degrees of freedom inside the cosmological horizon do not lead to diminishing cosmological 
constant. 

On the other hand, the degrees of freedom outside the cosmological horizon increase with 
time as more and more degrees go out of the horizon. Indeed we find that they may give 
rise to a desired effect. As they accumulate with cosmic evolution, the screening effect of 
the cosmological constant grows. Although it is suppressed by kA^, the suppression effect 
can be compensated by large logarithmic factors as arbitrary many degrees of freedom could 
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accumulate outside the horizon^ It is important to generalize our work to more generic 
situations including quantum effects of gravity. Since our investigation of the propagators 
is limited to the region well inside the horizon, we also need to understand the behavior of 
them around and beyond the cosmological horizon. 
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A Collision term evaluation 



In this appendix we explain the details of our calculation for the collision term. 



In the first step, using our integration formula (I4.10p . the on-shell collision term (14. 12^ is 
evaluated as 



Conlf] 

= + (1 + f{p))e-^^' 



X 



if 2 2 327rV JO 



-2f 



+ 



oo rpi+p 

dpi / dp2 

J\pi-p\ 



i{pi+P2-p) {pi+P2-pYr^ 

. 1 -2f -1 2 . 

i{pi -P2-p) (pi - p2 - py 

. 1 -2f -1 2 . 

^i{-pi+P2-p) {-Pi + P2 - pY 

. 1 -2f -1 2 . 

i{-Pi-P2-p) {-Pi - P2 - pY 

X {/(Pl)/(P2)-(l + /(Pl))(l + /(P2))} 



(A.l) 



■I- Such a mechanism is proposed by Tsamis and Woodard. See [H] and references therein. 
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- f{p) e 



+ipT 



+ { 



1 {-iXf 1 

1P 2 32^ JO 
1 +2f 



cxD rpi+p 
dpi / 

J |pi-p| 

-1 2 



+P2 -P) 



(Pl +P2 -PYt^ 



} 



+ {- 



+2f -1 2 ^ 



-P2-p) (pi -P2- pY 1"c 



+{ 



1 



-2r 



+ 



■} 



i{-pi+P2-p) (-Pl +P2 

. 1 +2f -1 2 . 

X {/(j>l)/(P2)-(l + /(j>l))(l + /(P2))} 



Here we have used the following relation. 

1 f d^pi d^p2 



2p J (27r)32pi (27r)32p2 

^ poo ppi+p 



P2 - Pj 



327rV 



dpi 



dp2. 



For the comparison with the off-sell part, we insert the identity factor as 

/ ^{27r)6{e-{±p,±p2)). 



(A.2) 



(A.3) 



In this way, we obtain the expression f l4.13p in the main text. The off-shell part is calculated 
just like the on-shell part. 



In the main text, we have introduced the collision terms with a finite energy resolution Ae 
following a standard procedure in massless field theories. They are given explicitly as follows 



C'Jf] ^ C^f] 



+ 



X 



IGirp'^H'^ 
de 



(A.4) 



27r 

Jp-Ae 
A2 



{e-py {e + pY r, 



2 



1 1 -1 

n — \2 - , X2 )~ / (1 + - ^) 



IGirp'^H'^ 
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"^^'ds 1 1 .-1 



p 



2 



2 



r dF 1 1 —1 /""^ 



~ de 1 1 .-1 



^0 

■X 



"^"'(7 ^ - r— / dp, f{p,)f{e - p,) 

+ 2 / / e+-^(- ^ - ^— ^)— / dp, f{p,){l + f{p, - s)) 

Jo 27r {e~py {e + py J e+e 



In the case of the thermal distribution function, the on-shell coUision term (]A.4p is evaluated 

as 



cLlf] 

A2 



,(l + f{p))e-^P^x 



IQupH'^ 

f+'^'.dE,, 1 1 ,-'(f , 1 1 ,-1-,,, „, 
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■X 



de 



+ 



27r 



1 + /(e))e-^ 
l + /(e))e- 



1 



1 



1 



-1 

I 



3a + G(£,p,/3)) 



de 



+ 



2tt 



p—Ae 



1 



{s + py^ 
1 



(A.6) 



G{e,p,f3) is defined in fl4.20p . We find that the hnear infra-red divergences at e = p are 
canceled, but the apparent logarithmic divergences remain. The situation here is analogous 
to QCD where the logarithmic divergences require the scale dependent modification of the 
parton distribution function. In our case, the IR singularity also leads to the modification 
of the particle distribution function as the final expression is shown in the main text f l4.22p . 



B Boltzmann equation in gcp^ theory 

In this appendix, we consider the Boltzmann equation in gip^ theory. Since this the- 
ory is classically stable, it is a good example for investigating quantum effects on the dS 
background. 

As in the main text, we evaluate the time integrations with the assumption \ {e ±p)rj| ^ 1. 
In Qip"^ theory, we need to retain higher order terms than f l4.10p to investigate the particle 
production effects in de Sitter space 



T 

^3 



1 —n —n(n + 1) 

A 1- 



n = 1,2,-- - . 



(B.i: 



We should note that (IB.ip can be evaluated exactly when n = 



oo 



dn e^("^P)"^ = e'("^P)"' X (B.2) 

i{€ ±p — iO) 

= e'^'^P^^^ X ( , ^ , + Tidie ± p)] . 

\i{e±p) J 
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The —iO prescription is necessary for the convergence at T3 = 



—00. 



In this appendix, we focus on the IR effects of the coUision term at £ — p = 0. Therefore we 
consider only 2 — > 2 processes. In these processes, the on-shell part of the coUision term is 
as follows 

Conif] = + (1 + m) e-^^' X / n {27rf6^'\p + Pi + P2 + Ps) 
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(B.3) 



The off-shell part of coUision term is as follows 
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(B.4) 
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In ( IB.Sp and (1B.4I) . only the leading term in 1/p|tc| expansion is shown for the energy 



conserving part containing 6{pi + P2 — P3 — p) ■ 

We note that the leading term in 1/p|tc| expansion is the same with the collision term in 
Minkowski space 

C[f] leading 

= y ^ / n J0^^ i^^rs^'^ (P + Pi + P2 + P3)5(P1 +P2-P3-P) (B.5) 

+ {/(Pl)/(P2)(1 + /(P3))(l + f{p)) - (1 + /(Pl))(l + f{P2))f{P3)f{p)} e-'^' 
- {f{Pl)f{P2){l + /(P3))(l + f{p)) - (1 + + f{P2))f{p3)f{p)} e+*^" 

This is because the leading term is conformally invariant. We thus obtain the identical result 
with [20] to the leading order in l/p|rc| expansion. 

In addition to the leading effect, we investigate the particle production effects due to energy 
non-conservation. Let us focus on the case that the initial distribution function is thermal. 
It solves the Boltzmann equation to the leading order as the following identity holds 

(1 + /(pi))(l + f{p2))fiP3)f{Pl+P2-P3) 
= f{Pl)f{P2){l + /(P3))(l + f{pi +P2- P3)). 
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Therefore the off-shell part is written as follows 

next leading 

- / n (2^)'^'='(p - p. + - P3) (B.T) 



X 



+ 3{(1 + + /(P2))/(P3) - /(Pl)/(P2)(1 + /(P3))} 

X (1 + /(pi + p2 - Ps)) e-'^'+P'-P^^^ 



+ P2 - P3 - P) P^ 
1 , 1 



1 


1 


1 




P2 


P3 


P 


c 


1 


1 


1 




— ^ 

P2 


P3 


p2 


c 


1 


1 




-4 


P2 


P3 


P 


' c 



«(Pl +P2 -P3 -P)^ 7>i 

+ 7 ^ ^ X — + 

(Pl +P2 -P3 -P)'' PI 

1 / 1 

+ 7 X — + 

(pi + P2 - P3 - pr Pi 

+ 3{(1 + + /(P2))/(P3) - /(Pl)/(P2)(1 + /(P3))} 

x/(pi+P2-P3) e+^(^i+P^-f^)^ 
. 1 1 -2f 

^ 1 '^n, -I- T)o — no — p2 ^3 

, 1 1 1 l,-2f 
X - + )^ 

Pi P2 P3 P Tc 



X 



1 


^(Pl +P2 - 
1 


-P3 


-P) 


?(Pl +P2 - 
1 


-P3 


_p)2 


(Pl +P2 - 
1 


P3 - 




(Pl +P2 - 


P3 - 


_p)3 



^ 1 1 1 1-2 

^" ' " " "^^ pf pI p2 ^3 

,1 1 1 1-4 . 

+ 7 X ( ^ )— [ 

(Pl +P2 -P3 -P)'^ Pl P2 P3 P 

Most of the IR divergences at pi + p2 — P3 — P = cancel out between Co„[/] and Co//[/]. 
This is because the total spectral weight is conserved due to unitarity. The remaining IR 
divergence comes from momentum dependence of the distribution function 

/(Pl+P2-P3)= /(P)H-/'(P)(P1+P2-P3-P) + --- ■ (B.8) 

As explained in the main text, this IR divergence leads to the change of the distribution 
function 

V ~ /' W / n (2-)'^'*(P + P: + p. + P3) (B.9) 
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X {(1 + + /(P2))/(P3) - /(Pl)/(P2)(1 + fiPs))} 




1. 2 1 



{Pl+P2-P3-Py Pi P2 P3 P rl .' 

Here again we may adopt the IR cut-off : |pi + P2 — Ps — p| ~ 1/|tc|- Tc dependence can be 
entirely absorbed into pliysical quantities Pi = piH\T^^ T = (3H\tc\. 

We may draw the following conclusion in this appendix. The leading order collision term 
is identical to that in Minkowski space. If we consider the higher order terms in 1/p|tc| 
expansion, the off-shell part is generated due to the particle production while the total 
spectral weight is preserved due to unitarity. We further find the non-trivial change of the 
distribution function due to IR divergences. These features in g(p^ theory are qualitatively 
identical to those in X(p^ theory. 

C Two and three point functions 

The contribution to the three point function from the degrees of freedom inside the 
cosmological horizon is evaluated as 
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It in fact does not possess IR singularity around £ — p ~ 0. 
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We subsequently evaluate the contributions from the degrees of freedom outside the horizon. 
We expand the original wave function in the power series of pr 

Ux) = -^^(1 + - i\{pTf ■ ■ ■ )e^^-\ (C.2) 

We also assume / = in these estimates. 

The expectation value of the three point function is 
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Here we have extracted only the leading term after rescaling —piT — > pi. 
The leading perturbative correction to the two point function is 
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-3 /p3X-3 



X (-r)-^ 

-3 



1 



l/max{pi} \ '1/ ••' — 

\ — 3 /n^ \ -3 



l/max{pi} ^2 



{(f) (f) -(f) (I) } 



l/max{pi} 

-3 



-'''1 J-l/mox{pi} (-"^2) 



{(f)-'(l)-'} 



X -T 



+ 



+ 



25 • 37r4 



l/max{pi} 

dp3 



dTi 



26 . 



/ C?P2 / ' 

J -EOT J —EOT 



;-log(-£oT))^ 



dcos6' 
1 P2P3 2 



l/max{pi} \ '2j 

1 1 



{\og{max{pi})y 



where p4 = -v/pi + p| + '^P2Pz cos 6'. 
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